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In this paper I intend to show that macroscopic entanglement is possible at high temperatures. I 
analyze multipartite entanglement produced by the n pairing mechanism which features strongly in 
the fermionic lattice models of high T c superconductivity. This problem is shown to be equivalent 
to calculating multipartite entanglement in totally symmetric states of qubits. I demonstrate that 
we can conclusively calculate the relative entropy of entanglement within any subset of qubits 
in an overall symmetric state. Three main results then follow. First, I show that the condition 
for superconductivity, namely the existence of the off diagonal long range order (ODLRO), is not 
dependent on two-site entanglement, but on just classical correlations as the sites become more and 
more distant. Secondly, the entanglement that does survive in the thermodynamical limit is the 
entanglement of the total lattice and, at half filling, it scales with the log of the number of sites. It 
is this entanglement that will exist at temperatures below the superconducting critical temperature, 
which can currently be as high as 160 Kelvin. Thirdly, I prove that a complete mixture of symmetric 
states does not contain any entanglement in the macroscopic limit. On the other hand, the same 
mixture of symmetric states possesses the same two qubit entanglement features as the pure states 
involved, in the sense that the mixing does not destroy entanglement for finite number of qubits, 
albeit it does decrease it. Maximal mixing of symmetric states also does not destroy ODLRO and 
classical correlations. I discuss various other inequalities between different entanglements as well as 
generalizations to the subsystems of any dimensionality (i.e. higher than spin half). 

I. INTRODUCTION 

Entanglement is currently one of the most researched phenomena in physics. Often shrouded in mystery, its basic 
premise is quite simple - entanglement is a correlation between distant particles that exists outside of any description 
offered by classical physics. Whilst this may at first glance seem an innocuous statement, in reality it is anything 
but. Predictions from the theory of entanglement have confounded some of the greatest minds in science. Einstein 
famously dubbed it spukhafte Fernwirkungen: "spooky action at a distance" . As we look deeper into the fabric of 
nature this "spooky" connection between particles is appearing everywhere, and its consequences are affecting the 
very (macroscopic) world that we experience. At an implementational level, using entanglement researchers have 
succeeded in teleporting information between two parties, designing cryptographic systems that cannot be broken 
and speeding up computations that would classically take a much longer time to execute [1]. Even though these 
applications have generated significant interest, I believe we have only scratched the "tip of the iceberg" in terms of 
what entanglement is, and indeed what we can do with it. 

Whilst entanglement is experimentally pretty much beyond dispute in microscopic systems - such as two photons 
or two atoms - many people find it difficult to accept that this phenomenon can exist and even have effects macro- 
scopically. Based on our everyday intuition we would, for example, find it very hard to believe that two cats or two 
human beings can be quantum entangled. Yet quantum physics does not tell us that there is any limitation to the 
existence of entanglement. It can, in principle and as far as we understand, be present in systems of any size and 
under many different external conditions. 

The usual argument against seeing macroscopic entanglement is that large systems have a large number of degrees 
of freedom interacting with the rest of the universe and it is this interaction that is responsible for destroying entan- 
glement. If we can exactly tell the state that a system is in, then this system cannot be entangled to any other system. 
In everyday life, objects exist at room (or comparable) temperatures so their overall state is quantum mechanically 
described by a very mixed state (this mixing due to temperature is, of course, also due to the interaction with a large 
"hot" environment). Mixing states that are entangled, in general, reduces entanglement and ultimately all entangle- 
ment vanishes if the temperature is high enough. The question then is how high is the highest temperature before we 
no longer see any entanglement? And how large can the body be so that entanglement is still present? Can we, for 
example, have macroscopic entanglement at the room temperature? 

Entanglement has recently been shown to affect macroscopic properties of solids, such as its magnetic susceptibility 
and heat capacity, but at a very low (critical) temperature [2] . This extraordinary result demonstrates that entangle- 
ment can have a significant effect in the macroscopic world. The basic reason for this dependence is simple. Magnetic 
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susceptibility is proportional to the correlation between nuclear spins in the solid. As we said before, entanglement 
offers a higher degree of correlation than anything allowed by classical physics and the corresponding quantum sus- 
ceptibility - which fully agrees with experimental results [2] - is higher than that predicted by using just classical 
correlations (for further theoretical support for this see my article in [3]). It is now very important to go beyond this 
low temperature regime and experimentally test entanglement at higher and higher temperatures. 

Thinking that high temperature entanglement is linked with (perhaps even responsible for) some other high tem- 
perature quantum phenomena, such as high temperature superconductivity, is tempting. After all, superconductivity 
is a manifestation of the existence of the off diagonal long range order (ODLRO) [4] which is a form of correlation 
that still persists in the thermodynamical (macroscopic) limit. However, it is not immediately obvious that this corre- 
lation contains any quantum entanglement. My main intention in this paper is to show that it does. This correlation 
contains multipartite entanglement between all electron pairs in the superconductor. To calculate this we need to be 
able to quantify entanglement exactly and be able to discriminate entanglement from any form of classical correlation. 

A great deal of effort has gone into theoretically understanding and quantifying entanglement [5]. There are a 
large number of different proposed measures; the different measures capture different aspects of entanglement. In 
this paper we will be interested in a measure that is based on the (asymptotic) distinguishability of entangled states 
from separable (disentangled) states known as the relative entropy of entanglement [6,7]. The main advantage of 
this measure is that it is easily defined for any number of systems of any dimensionality, which is not the case for 
entanglement of formation or distillation [5]. I have argued that a number of results in quantum information and 
computation follow from the relative entropy function [5] . 

There is, unfortunately, no closed form for the relative entropy of entanglement, but this measure can still be 
computed for a large class of relevant states such as the pure bipartite states, Werner states and many others [7]. 
Most recently, Wei et al [8] have succeeded in obtaining a formula for the relative entropy of entanglement for any 
number of totally symmetric pure states of n qubits using a very simple and elegant argument (some partial results 
have been obtained previously in this direction using different methods by Plenio and Vedral [9], but only for three 
qubit symmetric states) . I will use and extend these results further with the idea of applying them to a specific model 
of a superconductor. 

The purpose of this paper is to investigate possible links between high temperature entanglement and high tempera- 
ture superconductivity with the intention of showing that entanglement can persist at higher temperatures. I analyze 
a particular mechanism - the repairing of electrons due to Yang [10] - that was originally proposed to explain high 
temperature superconductivity. The chief difference between this pairing mechanism and the usual Bardeen, Cooper 
and Schrieffer (BCS) electron pairing [11] for (low temperature) superconductivity is that, in the former, electrons 
that are positioned at the same site are paired, while in the latter, electrons forming Cooper pairs are separated by a 
certain finite average distance (the so called coherence length, typically of the order of hundreds of nanometers). The 
physical reason behind electron pairing is also thought to be different in a high temperature superconductor, but I do 
not wish to enter into discussing these details here (see e.g. [12]). I will, however, look at the r\ model in a different 
way, using totally symmetric states, and this will make calculating entanglement easier. Wei et al. [8] have recently 
made very important steps in calculating the relative entropy of entanglement for symmetric state. I extend their 
approach to calculating the relative entropy of entanglement for mixed symmetric state arising from tracing over some 
qubits in pure states, and apply it to understanding various relations between entanglements of a subset of qubits 
and their relation to the total entanglement. I show that although two-site entanglement disappears as the distance 
between sites diverges (a conclusion also reached by Zanardi and Wang in a different way [13]), the total entangle- 
ment still persists in the thermodynamical limit. Furthermore, it scales logarithmically with the number of qubits. 
Therefore, it is this total entanglement that should be compared with ODLRO and not the two-site entanglement. 
While the two-site entanglement vanishes thermodynamically, two-site classical correlations are still present and so 
is the entanglement between two clusters of qubits (two cluster entanglement in 77 states has also been analysed by 
Fan [14]). I show that all aspects of my analysis can easily be generalised to higher than half spin systems. My hope 
is that this work - which is really just a first step in exploring high temperature entanglement - will be extended to 
different models with states other than symmetric and that this will allow us a much more complete understanding 
of entanglement and the role it plays in the macroscopic world. 

II. //-PAIRING IN SUPERCONDUCTIVITY 

The model that I describe now consists of a number of lattice sites, each of which can be occupied by fermions 
having spin up or spin down internal states. Let us introduce fermion creation and annihilation operators, c\ s and 
Ci, s respectively, where the subscript i refers to the zth lattice site and s refers for the value of the spin, | or J.. Since 
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fermions obey the Pauli exclusion principle, we can have at most two fermions attached to one and the same site. 
The c operators therefore satisfy the anticommutation relations: 

{ c i,s,c\ t } = 5ij5 Sit (1) 

and c's and c^'s anticommute as usual. (Some general features of fermionic entanglement - arising mainly from the 
Pauli exclusion principle - have been analysed in [13,15-17]). 

We only need assume that our model has the interaction which favors formation of Cooper pairs of fermions of 
opposite spin at each site [10]. The actual Hamiltonian is not relevant for my present purposes. It suffices to say that 
Yang originally considered the Hubbrad model for which the X] states are eigenstates (but none of them is a ground 
state [10]). A generalisation of the Hubbard model was presented in [18] and in a specific regime of this new model 
the r\ states do become lowest energy eigenstates (this is a fact that will become relevant when we talk about high 
temperature entanglement). Both these models have been used to simulate high-temperature superconductivity, since 
in high superconducting materials, the coherent length of each Cooper pair is on average much smaller than for a 
normal superconductor. 

Suppose, now, that there are n sites and suppose, further, that we introduce an operator r\ t that creates a coherent 
superposition of a Cooper pair in each of the lattice sites, 

n 

2=1 

The operator can be applied to the vacuum a number of times, each time creating a new coherent superposition. 
However, the number of applications, k, cannot exceed the number of sites, n, since we cannot have more than one 
pair per site due to the exclusion principle. I now introduce the following basis 

\k,n-k):=-^tf) k \0), (3) 



where the factor in front is just the necessary normalisation. Here, the vacuum state |0) is annihilated by all c operators, 
<k,s |0) = 0. Wc note in passing that the originally defined r\ operators can also have phase factors dependent on the 
location of the site on the lattice. We can have a set of operators like 

n 

All the states generated with any rjk from the vacuum will be shown to have the same amount of entanglement so 
that the extra phases will be ignored in the rest of the paper (i.e. we will only consider the k = states). 

We can think of the r\ states in the following way. Suppose that k = 2. Then this means that we will be creating 
two ry-pairs in total, but they cannot be created in the same lattice site. The state \2,n— 2) is therefore a symmetric 
superposition of all combinations of creating two pairs at two different sites. Let us, for the moment, use the label 
when the site is unoccupied and 1 when it is occupied. Then the state |2, n — 2) is 

|2,n-2) = — =(| 000 ...JTJ + -| JT, - 000)) (5) 

Y (2) n-2 2 2 n-2 

i.e. it is an equal superposition of states containing 2 states |1) and n — 2 states |0). These states, due to their high 
degree of symmetry, are much easier to handle than general arbitrary superpositions and we can compute entanglement 
for them between any number of sites. Note that in this description each site effectively holds one quantum bit, whose 
signifies that the site is empty and 1 signifies that the site is full. 

The main characteristic of 77 states is the existence of the long range off diagonal order (ODLRO), which implies 
its various superconducting features, such as the Meissncr effect and flux quantisation [19]. The ODLRO is defined 
by the off diagonal matrix elements of the two-site reduced density matrix being finite in the limit when the distance 
between the sites diverges. Namely, 
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where a is a constant (independent of n). I will show that although the existence of off diagonal matrix elements 
does not guarantee the existence of entanglement between the two sites, it does guarantee the existence of multi-site 
entanglement between all the sites. Note that here, by "correlations" I mean correlations between the number of 
electrons positioned at different sites i and j. Namely, we are looking at the probability of one site being occupied 
(empty) given that the other site is occupied (empty). This is different from spin-spin correlations, which would look 
at the occurrences of both electron spins being up or down, or one being up and the other being down [15]. 



III. GENERAL DESCRIPTION OF SYMMETRIC STATES 



The states I will analyze here will always be of the form 



|*(n,fc)) = \k,n- k) := 



(S\ 000 ... JA)) 

(fe) k n-k 



(7) 



where 5* is the total symmetrisation operator. We will also consider mixtures of these states, which become relevant 
when we talk about systems at finite temperatures. Symmetric states arise, for example, in the Dicke model in 
which n atoms simultaneously interact with a single mode of the electro-magnetic field [20]. They are, furthermore, 
very important as they happen to be eigenstates of many models in solid state physics, and, in particular, they are 
eigenstates of the Hubbard and related models supporting the r\ pairing mechanism. The analysis presented in this 
paper will be applicable to any of these systems and not just the r\ model. The r\ mechanism will be here significant 
because of its potential to support high temperature entanglement. 

I would now like to start to compute the entanglement between every pair of qubits (sites) in the above state 
\^(n,h)). A simpler task would be first to tell if and when every pair of qubits in a totally symmetric state is 
entangled. For this, we need only compute the reduced two-qubit density matrix which can be written as: 



a 12 (h) = a|00)(00| + &jll><llj + 2c|^ + )(^+| 



where = (|01) + |10))/\/2 and 



a = 



(fej) _ Hk - 1) 



(2) 

r k 2 ) 



n(n — 1) 

(n-h)(n-h-l) 



(1) ~ -("-1) 
(fc-i) = k(n-h) 

(2) n(n-l)- 



(8) 

(9) 
(10) 

(11) 



We can easily check that a + b + 2c = 1 and so the state is normalized. This density matrix is the same no matter 
how far the two sites are from each other, since the state is symmetric, and must therefore be identical for all qubits. 
We can easily test the Peres-Horodecki (partial transposition) condition [21] for separability of this state. Two states 
are entangled if and only if they are inseparable which leads to states a 12(h) being entangled if and only if 



a + b- v/(a-fe) 2 +4c 2 <0 



(12) 



which leads to 



(h- l)(n-h- 1) < h(n-k) . (13) 

This equation is satisfied for all n > 2 (two qubits or more) and 1 < k < n — 1. So, apart from the case when the total 
state is of the form |000..0) or 1 , there is always two-qubit entanglement present in symmetric states. Note, 
however, that in the limit of n and k becoming large - no matter what their ratio may be - the value of the left hand 
side approaches the value of the right hand side and entanglement thus disappears. This is a very interesting property 
of symmetric states and we will be able to quantify it exactly in the next section. 

An important point to make is that the two point correlation function used in the calculation of the ODLRO in eq. 
(6) is, in fact, just one of the sixteen numbers we need for the full two-site density matrix (the independent number 
of real parameters is actually fifteen, because of normalisation). In our simplified case of symmetric states in the 
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^-pairing model, this off diagonal clement is equal to c. However, for the density matrix we still need to know a and 
b, and these numbers clearly affect the amount of entanglement. Imagine, for example, the situation where a = b. 
Then the condition for entanglement is that a — c < 0, which does not hold if a > c and such a density matrix is 
certainly possible. So, the first lesson is that two-site entanglement is not the same as the existence of ODLRO, and 
therefore two-site entanglement is not relevant for superconductivity. This does not mean, of course, that there is 
no entanglement in the whole of the lattice. In the next section, I will calculate exactly this. We will determine the 
relative entropy of entanglement for all symmetric states and all their substates. I will be able to extend the method 
of Wei ct al [8] and analyze many relationships between various subsets of symmetric states, including the amount of 
entanglement in any subset of qubits (or sites). 



IV. RELATIVE ENTROPY OF ENTANGLEMENT FOR SYMMETRIC STATES 



The symmetric states are very convenient for studying various features of multipartite entanglement simply because, 
as we already indicated, we can compute exactly the relative entropy of entanglement for any reduced state including 
the total symmetric state for any n and k. It is expected that, because they possess a high degree of symmetry, 
they will also display a high degree of entanglement. It is precisely for this reason that they are suitable to allow the 
existence of entanglement at high temperatures. This will now be analyzed in detail. 

I first introduce the relative entropy of entanglement. The relative entropy of entanglement measures the distance 
between a state and the nearest disentangled (separable) state. If V is the set of all disentangled states (i.e. states of 
the form J2iPiPi ® 9%, ■ ■ ■ ® P l n i where Pi is any probability distribution), the measure of entanglement for a state a is 
then defined as 



£(a):=mm S{a\\p) , (14) 

p62? 

where S(o"||/o) = tr(a log a — alogp) is the relative entropy between the two density matrices p and a. In order to 
compute this measure for any state a we need to be able to find its closest disentangled state p. Finding this closest 
state is, in general, still an open problem, however, it has recently been solved for pure symmetric states by Wei et al 
[8]. 

Wei et al showed that a convenient and intuitive way of writing the closest disentangled state to the symmetric 
state \k, n — k) is [8]: 



i-^ 27r ^|^«)(^»|, (15) 



2 :v ' 

where 

= (v/fcT^O) + y/(n - k)/ne l ^\l))® n (16) 

is the tensor product of n states each of which is a superposition of the states |0) and |1) with probabilities kjn and 
1 — k/n respectively. This p was proved to achieve the minimum of the relative entropy by showing that it saturates an 
independently obtained lower bound. The relative entropy of entanglement of the total state is now easily computed. 
Since a — |fc,n— fc)(fc,n— k\ is a pure state, traloga = and we only need to compute — (k,n— k\\ogp\k,n — fc), 
which is equal to 

(77 \ 77 77 

J + fcIog- + (n -*)Iog— . (17) 

Note that entanglement is largest when n = 2k as is intuitively expected (i.e. the largest number of terms is then 
present in the expansion of the state in terms of the computational basis states). Then, for large n, it can be seen 
that the amount of entanglement grows as 

£(|n/2,n/2»w^(logn + 2) (18) 

and so (in the leading order) entanglement grows logarithmically with the number of qubits in the state. To obtain 
this formula I have used Sterling's approximation for the factorial 
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n\ w 2.507n™ +1 / 2 e"™ . (19) 

Most results in this paper will asymptotically have the form a log n + (3 where a > and (3 are constants that will 
usually be omitted as we only care about the general form of the behaviour. 

I now return to the question of different phases introduced between different elements of the superposition in the 
symmetric states. Let us consider states of the form 

|1, n - 1, 6) = |00..1) + e* e |00..10) + e^^lCO) , (20) 

where we have k = 1 ones and n — 1 zeroes and 8 is any phase. The simplest way of seeing that entanglement does 
not depend on the phase 9 is to define a new basis at the mth site as |0) = |0), |1) = exp{(m — 1)0}|1). This way the 
phases have been absorbed by into the basis states and the resulting state is, in the tilde basis, 

|1, n-l,6} = |o6..i) + |66..l6) + |10..6) . (21) 

The amount of entanglement must therefore be independent of any phase difference of the above type and this is, of 
course, true for symmetric states with any number of zeroes and ones. All considerations from this point onwards will 
therefore immediately apply to all these states will different phases. 

We can also compute the two-site relative entropy of entanglement exactly. The closest disentangled state is in this 
case the same as in eq. (15) with n = 2. In the computational basis we have 

p = |00><00|+ (^) 2 |11)(H|+ {^T 3 ) • (22) 

That this is a minimum can be seen from the fact that the relative entropy of the state of two qubits is: 

S{a\\p) = -S(a)-^ + \\ogp\^+) - (00|logp|00) - (ll|logp|ll) (23) 
> -S(a) - log<V + |/C# + ) - log(00|p|00> - log(ll|p|ll) , (24) 

the inequality following from concavity of the log function. Suppose now that p's only non-zero elements are poo — 
(00|p|00), pn = (ll|p|ll) and p ++ = {ip + \p\ip + ). Given that it has to be separable, meaning that 2 % / / o oPii > P++ 
(which follows from the Peres-Horodecki criterion), and that, at the same time, it has to be closest to a, we can 
conclude that poo = k/n. The other entries of p then follow. 

To prove that p is the minimum in a rigorous fashion, we need to show that any variation of the type (1 — x)p + xu> 
where u> is any separable state leads to a higher relative entropy (a method similar to [7] ) . Since relative entropy is a 
convex function, this means that 

-^S(a\\(l-x)p + xcj) >0. (25) 
dx 

In fact, since relative entropy is convex in the second argument it is enough to assume that w is just a product state. 

For a > 0, loga = J °° ^ffj^i, and thus, for any positive operator A, log A = J °° ^nTTrl 7 ' f^^) = 
S(a\\(l- x)p + xuj). Then 

^(0, W ) = - lim Tr (^g((l-^ + ^)-logP) 

OX x^Q I X 

= Tr|(a^ (p + t )-\p-u){p + t)- x dt)} 

= i-\ r TY{o{p + ty 1 bj{p + ty 1 )dt 

Jo 

r-OQ 

= 1- Tr((p + t)- 1 cr(p + t)- 1 uj)dt . (26) 
Jo 

For our minimal guess p in eq. (22) we can then write 

g(o,c) - 1 = -ii{u> J™ {p + ty^ip + t)- 1 ^ 

n k — 1 ,„„. n n — k— 1, . . 

+ , (27) 

n — 1 



G 



where we have used the fact that J^°(p + t) 2 dt = p 1 . Since the expression in the previous equation is always less 
than or equal to a unity if u = \ap)(a(3\ (i.e. a product state), it follows that 



< 1 • (28) 



Thus it also follows that §j(0, \a/3)(a/3\) > 0. But any separable state can be written in the form p = ^\ r^a 1 (3 % ) (a* /3 Z 
and so 



g(0,p) = 5>§£(°' > . (29) 

i 

And this confirms that p is the minimum since the gradient is positive. 
Therefore the relative entropy of entanglement between any two sites is: 

E\2 = a log a — b log b — 2c log 2c 

, f k V fn-k\ 2 n , (2k(n-k) 

-a log - -6 log -2clog v ; 

\n/ \ n J \ n z 

, / n \ fc(fc-l), (n-fc)(n-fc- 1) , /n-fc-l\ 

We see that when n, fc, n— fc — > oo, then £12 — > as it should be from our discussion of the separability criterion. This 
can be thought of as one way of recovering the "quantum to classical" correspondence in the limit of large number of 
systems present in the state: locally between any two sites, entanglement does vanish, although globally, and as will 
be seen in more detail, entanglement still persists. 

Entanglement of any number of qubits, I < k, can also be calculated using the same method . The state after we 
trace out all but I qubits is given by 

*i=YJ, 1 (31) 



^V-iJ (I) 

The closest disentangled state is given by 



« = E( U) (V^lM-w-ii, 02) 



as can be shown by the above method. The relative entropy of entanglement is now given by 



This is a very interesting quantity as it allows us to speak about entanglement involving any number of qubits. What 
do we expect from it? We expect that entanglement grows exponentially with /, for a fixed total number of qubits, n. 
This can be confirmed using the Sterling formula. Note that entanglement grows at this rate even though the states 
we are talking about are mixed, since n — I qubits have been traced out. Another way of seeing why entanglement 
grows exponentially with the number of qubits included for a total fixed number of qubits, is to look at the opposite 
regime. For any finite fixed /, we should have that in the large n,k,n — k limit the amount of entanglement between 
/ tends to zero. This decrease with larger and larger n happens at an exponential rate. 



V. CLASSICAL VERSUS QUANTUM CORRELATIONS 

In this section I would like to investigate the relationship between classical and quantum correlations for symmetric 
states, and both in relation to the already introduced concept of ODLRO. First of all, it is clear that in the limit of 
n — > 00 all bipartite (or two-site) entanglement disappears (this was seen both from the Peres-Horodecki criterion and 
from the direct computation of the relative entropy). In spite of this, the ODLRO still exists and the two quantities 
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are therefore not related. In other words, two-site entanglement is not relevant for superconductivity. However the 
main point of this section is that the two-site classical correlations still survive in the limit of n — > oo. In order to 
show this, let us, first of all, define bipartite classical correlations. 

A quantum state can have zero amount of entanglement, but still have non-zero classical correlations. An example 
is the state 1 00} (00 1 + |11)(11|. Classical correlations between systems A and B in the state gab can be defined as 
[22] 

Ca{<?ab) ■= maxSVs) - VVS(ctb) = max ^PiS^BW^B) , (34) 

A'.Ai . A^.Ai^ 

i 1 % i 1 % 

where g b — trA<J AB , g ab = A^abA], and X^i = 1 is the most general measurement on system A. The same 

can be defined with the most general measurement performed on B, so that we obtain 

C b ((tab) ■= max S(a A )-y^PiS(a A ) = max V Pi S(g a \\g a ) ■ (35) 

B!Bi B'Bi- 

The physical motivation behind the above definition is the following: classical correlations between A and B tell us 
how much information we can obtain about A (B) by performing measurements in B (A). It is the (maximum) 
difference between the entropy of A (B) before and after the measurement on B (A) is performed. There is some 
evidence that Ca = Cb [22], but this equality will not be relevant here. 

Now, applying this measure of classical correlations to the two-site reduced density matrix from the overall sym- 
metric state, pi2, we obtain, 

C = -aloga- b\ogb- clogc+ i((a + c/2) log(a + c/2) + (b + c/2) log(6 + c/2)) 

= (r - 2r 2 ) log r + ((1 - r) - 2(1 - rf) log(l - r) - 2r(l - r) log 2r(l - r) (36) 

where r = k/n is the fraction of ones in the state (the so called filling factor in any "Cooper pair" lattice model, 
including the r\ model). We now see that at half filling - when ODLRO is maximal - the classical two-site correlations 
also survive asymptotically since Ca = Cb = 0.5. Therefore, all the correlations between any two sites are here due 
to classical correlations. 

Note, incidentally, that we cannot have the situation in which entanglement exists between two parties, while at the 
same time classical correlations vanish. Quantum correlations presuppose the existence of classical correlations. This, 
of course, relies on the fact that entanglement is defined in a reasonable way, namely that when we talk about two-site 
entanglement we must trace the other sites out. We are not allowed to perform measurements on other sites and 
condition the remaining entanglement on them. Measurements that generate entanglement are, first of all, unrealistic 
for a macroscopic object which thermalizes very quickly. Even if we were to allow such measurements, then the state 
after them will still have classical correlations of at least the same magnitude as entanglement. So, it cannot be that 
entanglement is important for the issues of superconductivity, phase transitions, condensation, etc., and that classical 
correlations are not. 

As an example, let us take the "maxmum singlet fraction" in the two-site density matrix g\i as our definition of 
entanglement. This is the maximaum fraction of a maximally entangled state in the state G\2, which is in this case 
equal to c, and this is the same as ODLRO. So, if the maximum singlet fraction is used to measure entanglement, then 
entanglement also persists in the thermodynamical limit. In fact, as will be shown later, this measure also survives 
when we mix symmetric states, because it is a linear measure. The maximum singlet fraction, however, is not a 
realistic measure of entanglement as it is not easily accessible experimentally, which is why we do not use it in this 
paper. 

In order to make our analysis more complete we also show how to calculate mutual information [5] for symmetric 
states. This quantity tells us about the total (quantum plus classical) correlations in a give state. Mutual information 
is equal to the relative entropy between the state itself and the product of individual qubit density matrices, obtained 
by tracing out all the other qubits. This product state is easily written down to be: 

_| )(0| + ^— |l)(l|j . (37) 

The mutual information is now given by 

I{\k,n-k))=n[ — log log , 38 

\ n n n n J 
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and this is basically just the sum of individual qubit entropies. Since the qubit entropy (the quantity in brackets in 
the above equation) is a finite quantity for a given ratio r = k/n, the total mutual information grows linearly with 
the number of qubits n. Furthermore, since entanglement grows as log n, we conclude that classical correlations grow 
roughly as n — logn (for this conclusion to be exact, classical and quantum correlations as defined here would have 
to add up to mutual information; while this is true for some states [22], it is certainly not true in general). 

The fact that classical correlations and mutual information survive the thermodynamical limit does not imply that 
there in no meaning left for entanglement when it comes to superconductivity and ODLRO. Only now, we must talk 
cither about the bipartite entanglement between two clusters of sites (to be computed in the next section) or the 
multipartite entanglement between all sites. Since the overall state across all sites is pure in our considerations so 
far, this means that two-site non-vanishing classical correlations (or equivalcntly ODLRO) must imply entanglement 
between two clusters, each of which contains one of the sites and such that the union of the two clusters is the whole 
lattice. This simply must be the case, since, otherwise, if the clusters were not entangled, the total state would be 
a product of the states of individual clusters, and this means that even classical correlations would be zero, which 
is a contradiction. Furthermore, the fact that any two such clusters are entangled, must mean that the multipartite 
entanglement also exists, since this entanglement is by definition larger than any bipartite entanglement (as, for 
multipartite entanglement, we are looking for the closest separable state over all sites, rather than just over the two 
clusters). I now quantify these various relations a bit more precisely. 



VI. VARIOUS OTHER RELATIONS BETWEEN ENTANGLEMENTS 



In this section I will discuss some other results that can be derived from our knowledge of symmetric states so 
far. Some of the results will not necessarily be directly relevant for the main theme of the paper - high temperature 
entanglement - but this section is a natural place to present them. The discussion about high temperature entan- 
glement in the rest of the paper can be understood without reading this section. I will fully return to the main 
topic in the next section. The first important question to be addressed here is the following. Suppose we look at the 
entanglement between one qubits and the rest n — 1 qubits in total and individually. We would expect that the total 
one-versus-rest entanglement is larger then the individual sum of the two-qubit entanglements. The logic behind this 
conclusion is that by looking at entanglements individually we always lose something from the total entanglement, 
i.e. the operation of tracing reduces entanglement. This translates into the following inequality: 

(n - 1)^12 < £l:(2,3..n) , (39) 

where 

„ k, f k\ n — k, (n—k\ 

£l:(2.3..n) = — fog fog 40 

v 1 n \n J n \ n J 

is basically the same as the entropy of every qubit in the symmetric state. We can prove this inequality by noting 
that it holds for k = 1 and 2k = n (the extreme points) , and because of the monotonicity and continuity of both sides 
it has to hold in general. 

The aforementioned inequality has a very important implication which shows that the bipartite entanglement in 
the symmetric state is always bounded from above by 

E l2 <^^<^-~ 1 - (41) 
(n — 1) n — 1 n 

the second inequality following from the fact that the entanglement between one qubit and the rest is equal to the 
entropy of that qubit and that can at most be log 2 = 1. Therefore, while the total entanglement of the symmetric 
state increases with log n, the two qubit entanglement decreases as 1 /n. Here we see most directly how it is possible to 
have the emergence of (only) classical correlations between constituents even though globally entanglement increases. 
There are many other open questions related to this one. We can repeat the same calculation for any fixed number of 
qubits. We can check if a cluster of qubits is, for example, more entangled to another cluster of qubits in total or if 
we add all the entanglements between their individual elements. Some of these may not be easy questions to answer 
in general. 

I would now like to calculate the entanglement between I qubits and the remaining n — I qubits. Since the whole 
state that we are now examining is pure, the relative entropy of entanglement is given by the entropy of the I qubits: 
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What are the properties of this expression when we take the various asymptotic limits? How is this quantity related to 
other entanglements calculated here? We expect that for the half filling, n/k = 2, and ra, I — > oo, the entropy becomes 
log/, since we basically have a maximal mixture in the symmetric subspace of I qubits. This can be confirmed by a 
simple application of the Sterling approximation formula used before. The result is in agreement with the fact that 
total entanglement grows at the rate of the log of the number of qubits, since two cluster entanglement is a lower 
bound for the total entanglement in the state between all the qubits. 

The last question I address is the relationship between the lower and higher order entanglement in the symmetric 
states. More precisely, the question is: if we add all the entanglements up to and including m qubits, is this quantity 
larger or smaller than the amount of entanglement of to + 1 qubits? Mathematically, this translates into the following 
two possible inequalities: 

m m 

Ei < E m+1 or E i > E m+i (43) 

i=l »=1 

where E m is given in cq. (33). We already know that for n = 3 and k = 1, 2, and I = 3 we have the equality in the 
above, namely E 3 = E\ + E 2 [9]. From this result alone it is not clear which way to expect the inequality to be. 
Numerical examples show us that, in fact, both results are possible. If we check the inequality for n — 100, k = 50 
and I = 4 for example, than the left hand side is smaller than the right hand side and the first inequality holds. For 
n = 100, k = 50 and I — 30, on the other hand, the left hand side is larger than the right hand side and the second 
inequality is satisfied. It is an interesting and open question to investigate the point of the cross-over when the two 
sides become equal to each other. 



VII. THERMAL ENTANGLEMENT AND SUPERCONDUCTIVITY 



There is a critical temperature beyond which any superconductor becomes a normal conductor. The basic idea 
behind computing this temperature according to BCS is the following. At a very low temperature, only the ground 
state of the system is populated and for a superconductor this state involves a collection of Cooper pairs with 
different momenta values around the Fermi surface. This state can be, somewhat loosely, thought of as a Cooper pair 
condensate, and it is this condensation that is the key to superconductivity. It took initially a long time to understand 
how the pairs are formed, since electrons repel each other and therefore should not be bound together. The attraction 
is provided by electrons interacting with the positive ions left in the lattice. We can think of one electron moving and 
dragging along the lattice, which then pulls other electrons thereby providing the necessary attraction [11]. When the 
temperature starts to increase, the Cooper pairs start to break up, leading to the transition to the normal conductor. 
What this "breaking up" means is that higher than ground states start to get populated by electrons, and these are 
states where an electron is created with say momentum k and spin up, but no electron is created in the — k momentum 
state. From the BCS analysis this critical temperature can be calculated to be of the form [11] 

T^^e-VA (44) 

where Tiu is the energy shell around the Fermi surface which is engaged in formation of Cooper pairs, k is the 
Boltzmann constant and A is a parameter equal to the product of the electron density at the Fermi surface N(0) and 
the effective electronic attractive coupling, V. The critical temperature formula is valid in the weak coupling regime 
where A = N(0)V « 1. 

The formula for the critical temperature is usually used for other mechanisms of electron pairing, and not just 
coupling via the phonon lattice modes as in the BCS model [11]. Importantly for us, the formula also features in 
models for explaining and designing high T c superconducting materials. If the attraction, say between an electron 
and a hole, is of the order of Coulomb forces, Uu> w leVm, and for the weak coupling of, say, A = 0.2, the critical 
temperature we obtain is 100K . So if the material is below this temperature, it is then superconducting. Anything 
above 70 — 90K is considered to be high temperature superconductivity, since it can be achieved by cooling with 
liquid nitrogen (which is a standard and easy method of cooling). What seems to be the mechanism behind high 
temperature superconductivity, is the fact that the energy gap between the ground superconducting, electron-pair 
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state, and the excited states is large enough not to be easily excited as the temperature increases well beyond zero 
temperature. The exact way in which this is achieved is still an open question. In the models mentioned here the 
ground state is one of the symmetric states from the previous sections. Therefore, we can conclude that as long as we 
have high temperature superconductivity, the total state should also be macroscopically entangled. Superconduction 
and hence entanglement can currently exist at temperatures of about 160 Kelvin. 

I would now like to explicitly calculate and show how entanglement disappears as the temperature increases for 
any model having the r\ pairing state as the ground state. For this, we need to be able to describe other states that 
would be mixed in with the symmetric r\ states as the temperature increases. They, of course, depend on the actual 
Hamiltonian. For instance, in the Hubbard model in [10], states of the type 

^l°)=E c ll c ! +a ,Tl°) ( 45 ) 

i 

are important; here we create a spin singlet state but at sites separated by the distance a^O. If we have 2k electrons 
in total, then 2k — 2 would be paired in the lowest energy state, and the remaining two electrons would not be. This 
would give us the state of the form 

10 ^-^lO). (46) 

Note that this state is a symmetric combination of states which have k — 1 electron pairs distributed among n sites and 
the last electron pair is in two different sites separated by the distance a. These two sites are different from the other 
k — 2 sites due to Pauli's exclusion principle. Even higher states are obtained by having two electron pairs existing 
outside of the symmetric state and so on. The exact form of these, as noted before, depends on the exact form of 
the Hamiltonian. Even simple Hamiltonians are frequently very difficult to diagonalize and their eigenstates are still 
by and large unknown. Given this, it may be difficult to calculate the exact amount of entanglement when, at finite 
temperature, the ground state is mixed with higher energy states. I will, therefore, make a simplifying assumption 
that, if the ground state is \k, n — k), the higher energy states can be written as \k — 1, n — k + 1), \k — 2, n — k + 2) 
and so on. All these will in fact be assumed to be symmetric and I will ignore the extra unpaired electrons as far as 
entanglement is concerned (they will only contribute to the eigenvalue of energy as it were) . 

This assumption leads us to consider mixtures of symmetric states. The symmetric states will be mixed with 
probabilities in accordance to Boltzmann's exponential law, or the Fermi-Dirac law if we talk about ry pairs. Which 
distribution we use will be immaterial for our argument. The total state, or, is 

n 

a T = $>*!*(*» »)><*(*> »)l (47) 

fc=0 

where, in the case of r\ pairs, the probabilities are 

Pi = £ eJt + I ( 48 ) 
where pi is the probability of occupying the ith energy level. The reduced two-site state can be calculated to be 

n 

V12 = ^Pk<y\2{k) . (49) 

fe=0 

The condition for inseparability now becomes 

^2 PkP ik(n - l){(n - k)l - (k - l)(n - I - 1)} > . (50) 
k.i 

We see that the thermal averaging is in a sense inconsequential for the existence of entanglement as the factors pkPi 
are probabilities and are always non-negative. For inequality to hold (i.e. to have non-zero bipartite entanglement 
present) we need that 1 < k,l < n — 1. This is the same condition as before when the total state was pure. Thus, 
surprisingly, the condition for inseparability is completely independent of temperature (although, two-site states do 
become separable in the macroscopic limit even at zero temperature, as noted before). 

We now look at the entanglement of the symmetric mixed state as a whole. Can we still calculate the relative 
entropy of entanglement? This is in general very difficult to do for multiparty mixed states, and some partial methods 
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for upper bounds have only been presented recently [3]. We conjecture that the closest disentangled state is now 
presumably the thermal average of the closest disentangled states for individual fc's (this, I believe, is the same as 
the conjecture in [8], for which Wei et al have offered a great deal of "circumstantial evidence"; for example, closest 
separable states have to possess the same symmetry as the entangled states for which they minimise the relative 
entropy [23]). I believe that this bound is exact and that this can be proven using methods for calculating two site 
entanglement, but I have not been able to show this yet. Even if this is not true, my method at least gives us a 
very good upper bound which is sufficient to show how total entanglement vanishes as T becomes high. The relative 
entropy of entanglement between these two states is given by (the right hand side of the inequality) 

E{a T ) < J> lcg Pfe - 5><*(M)| l og Y^PWl !*(*,»)> (5!) 

k k \ I ) 

where pi is the closest disentangled state to the pure symmetric state containing I ones and n — I zeroes. We have 
already seen that 

« = E (!) Ml,n))(n,n)\ , (52) 



i=l 

so that 



l\ ( fc V 4 I n — fc" 



E(a T )< ^Pk logpfe - 

k 

- »)|l°g{l>£Q (yj [^) Jl*(f,»))<*(l,n)ll*(*,»)> 



i=l 



The interesting conclusion here is the following. Suppose that we arc at a high temperature and that all symmetric 
states are equally likely, meaning that pt = l/(n + 1) for all values fc (basically, our state is an equal mixture of all 
symmetric states). This, of course, in an approximation to the true density matrix, but it becomes more and more 
accurate with the increase of temperature and it ceases to be so when states other than symmetric become mixed in. 
The (upper bound to) entanglement is then given by 

1 , \ /fc\ f k\ k 1 ( n — k 



%—\ 

The fraction inside the log tends to n 2 as n becomes large, so that entanglement scales as logn/n. This is to be 
expected as entanglement grows as logn with n, but the mixedness grows linearly with the number of state involved, 
n + 1. Therefore, in the thcrmodynamical limit, the overall mixed state entanglement also disappears. This has 
to eventually happen, of course, if we believe that entanglement is intimately linked with superconductivity and 
superconductivity also vanishes at sufficiently high temperatures. 

One kind of entanglement that we can say survives the thcrmodynamical high temperature limit is the average of 
entanglements of individual symmetric states. This average entanglement is given by 

E avr = J2PkE(\Kn - fc)) = \ ^> fc k>g Hn ~ k) . (55) 

k k 

Note that if all probabilities go as 1/n - i.e. the symmetric state is maximally mixed, then the entanglement scales as 
log n (the same as pure state at half filling) . This is expected, as there are n+1 states and each one has entanglement 
proportional to logn and, so, on average, entanglement also goes as logn. However, this average entanglement, as 
we argued before is not a good measure as it requires us to be able to address the symmetric states individually and 
discriminate them from each other. This is not just difficult in practice, but is in fact frequently even impossible in 
principle. 

It is interesting to note that the ODLRO does survive the mixing of symmetric states. Even when we have an equal 
mixture of all symmetric states the average ODLRO is given by 
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1 y^fc(n-fc) _1 12n+l 1 
n+ 1 2-" n n ~2~6 n ^6 ^ ' 

where the arrow indicates the convergence when n is large. Of course, at sufficiently high temperatures the system 
will leave the subspace of symmetric states and other states will also start to contribute. This eventually does lead 
to vanishing of ODLRO, but the total entanglement and ODLRO may still disappear at different temperatures. To 
calculate this exactly, we would need a much more detailed model and a more extensive and careful calculation which 
lie outside of the scope of the present paper. (Note: the same conclusions hold for the maximum singlet fraction in 
the two-site density matrix which also survives the mixing in the thermodynamical limit; this is, unfortunately and 
as pointed out before, not a suitable measure of entanglement in our setting). 

I conclude by showing that total correlations - quantum and classical - as quantified by the mutual information [5] 
can also easily be calculated for thermal mixtures of symmetric states. Let us assume again that the symmetric states 
are maximally mixed and each appears with the probability l/(n+ 1). Then the mutual information is given by 

I= — -V log log -logn+1 . 57 

n + 1 V n n n n I 

k x ' 

For large n this expression reduces to 

I — > n — \ogn . (58) 

Since wc know that thermal entanglement disappears in this limit, it is natural that the mutual information is equal 
to the classical correlations and this then coincides with the conclusion following eq. (38). 



VIII. D-DIMENSIONAL SYMMETRIC STATES 

Extensions of all the considerations in this paper to D-dimensions are seen to be very straightforward (similar 
generalizations to higher dimension symmetric states were also considered in [8]). We should actually be able to 
reproduce all the above results in the generalized form, such that instead of qubits we have qutrits, and so on. The 
generic symmetric state would now be written as: 

\ni,n 2 ,..n d ) (59) 

and it would be a totally symmetrized state of n\ states |1), ni states |2) and so on (it is also realistic to assume that 
the total number of particles is conserved). This could, for example, represent higher spin fermions which can occupy 
different lattice sites as in the rest of the paper. The closest state to the one in eq. (59) in terms of the relative 
entropy is a mixture of the states of the type 

(v^V|l) + e*yfc/N\2) + ...e^-WV^VId))®" , (60) 

with the phase <j> completely randomized as before. Knowledge of this closest state allows us to compute the relative 
entropy of entanglement of any number of subsystems of this system. All other results follow in exactly the same way. 
Nothing fundamental is changed in higher dimensions which is why I will not say anything more on this topic. 



IX. DISCUSSION AND CONCLUSIONS 



In this paper I have analyzed the rj pairing mechanism which leads to eigenstates of the Hubbard and similar 
models used in explaining high temperature superconductivity. I have shown that they correspond to multi qubit 
symmetric states, where the qubit is made up of an empty and a full site (two-electron spin singlet state). I have also 
shown how to calculate entanglement and classical correlations for such states. For pure states, entanglement of the 
total state increases at the rate logn with the number of qubits n, while two-site entanglement vanishes at the rate 
1/n. The two-site classical correlations, on the other hand, persist in the thermodynamical limit. So, the ODLRO 
can be associated for pure states with total entanglement or two-site classical correlations, but not with the two-site 
entanglement. I have also demonstrated that the total entanglement for maximally mixed symmetric states disappears 
at the rate (logn)/n. Various mutual information measures, which quantify the total amount of correlations in a given 
state, are also computed and shown to be consistent with the calculations of classical and quantum correlations. 
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There are many interesting issues raised by this work. Even if a consensus is reached on the correct model for high 
T c superconductivity, and this is shown to contain multipartite electron entanglement - which we have argued for in 
this paper - we are still left with the question of being able to extract and use this entanglement. At present there are 
no methods of extraction. Perhaps we can somehow extract electrons from the superconductor and then use them for 
quantum tcleportation or other forms of quantum information processing. 

It is presently believed that in order to perform a reliable and scalable quantum computation we may need to be 
at very low temperatures, but the existence of high temperature macroscopic entanglement may just challenge this 
dogma. Be that as it may, I believe that the argument in favor of the existence of high temperature entanglement 
does show that entanglement may be much more ubiquitous than is presently thought. This may force us to push the 
boundary between the classical and the quantum world towards taking seriously the concept that quantum mechanics 
is indeed universal and should be applied at all levels of complexity, independently of the number or, indeed, nature 
of particles involved. 
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